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Abstract 

The ever increasing adoption of mobile devices with limited energy storage capacity, on 
the one hand, and more awareness of the environmental impact of massive data centres and 
server pools, on the other hand, have both led to an increased interest in energy management 
algorithms. 

The main contribution of this paper is to present several new constant factor approximation 
algorithms for energy aware scheduling problems where the objective is to minimize weighted 
completion time plus the cost of the energy consumed, in the one machine non-preemptive 
setting, while allowing release dates and deadlines. Unlike previous known algorithms these new 
algorithms can handle general job-dependent energy cost functions, extending the application 
of these algorithms to settings outside the typical CPU-energy one. These new settings include 
problems where in addition, or instead, of energy costs we also have maintenance costs, wear and 
tear, replacement costs, etc., which in general depend on the speed at which the machine runs 
but also depend on the types of jobs processed. Our algorithms also extend to approximating 
weighted tardiness plus energy cost, an inherently more difficult problem that has not been 
addressed in the literature. 
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1 Introduction 

Managing energy consumption is a problem of critical interest throughout the world and throughout 
various industries. Computing devices use a large amount of energy, both in individual devices such 
as laptops and PDAs and also in large industrial uses such as datacenters. For example, Google 
states that the servers in its datacenter, which are much more efficient than the average industry 
server, consume IkJ per query on average [1]. In January 2011, just in the US, there were an average 
more than 400 million queries per day [2], and thus the total amount of energy consumed was 44.5 
million kWh, equivalent to more than 4, 000 average US households |3j. Furthermore, CPUs account 
for 50-60% of a typical computer's energy consumption [3], making CPU energy management very 
important. When scheduling on such devices, it is important not only to consider the relevant 
quality of service (QoS) metrics such as makespan or weighted completion time, but also to take 
energy consumption into account. Most modern CPUs can be run at multiple speeds; the lower 
the speed, the less energy used, and the relationship is device-dependent, but typically superlinear. 
The technique of scheduling while controlling the speed of the processor is known as speed scaling. 

Starting with the work of Yao, Demers, and Shenker [25], there has by now been tens of papers 
studying scheduling problems in which energy consumption is taken into account. (See, for example, 
the surveys by Irani and Pruhs |19] and that by Albers jl]). There are three main settings for energy 
aware scheduling problem: optimizing a QoS metric with an energy budget \22\ 123] . minimizing 
energy subject to a QoS constraint [71 [TOl [lU [25], or optimizing some convex combination of a 
scheduling objective and energy consumption [5l [6l [9l [12] . Underlying the latter setting, which is in 
the one we will focus in this work, is an assumption that both energy and time can be (implicitly) 
converted into a common unit, such as dollars. 

1.1 Our results 

In this paper we consider two commonly studied scheduling metrics, weighted completion time 
and weighted tardiness, that have not received attention in the energy aware scheduling literature. 
Given a schedule in which job i with weight Wi, release time r^, and deadline di is completed 
at time Cj, the total weighted completion time is ^^WiCi. The tardiness of a job is zero if it 
is completed before its deadline and otherwise equal to the amount by which it misses, that is, 
Ti = max{0, Cj — di\ and total weighted tardiness is ^ WiTi. For both these metrics, we consider 
the non-preemptive, off-line problem on one machine, and allow arbitrary precedence constraints. 
For the weighted completion time we allow arbitrary release dates as well. We consider a metric 
that is a convex combination of our scheduling metric and energy cost. We are not aware of any 
previous work on energy aware scheduling algorithms for these metrics. There is a rich literature 
on minimizing weighted completion time in the absence of energy concerns (e.g. |20 1 [211 IM] ). but 
we are aware of only one result about weighted tardiness in the absence of energy concerns in the 
speed scaling/resource augmentation literature [8], where a 2-machine, 24-speed 4-approximation 
algorithm is presented. Weighted tardiness, in particular, is difficult to analyze because, in contrast 
to most scheduling objectives, it is a non-linear function of completion time. 

In our work we consider a more general model of energy cost than has previously been used. 
The most common energy model assumes that the rate at which power is consumed is a polynomial 
function of speed of the form P{s) = s^ for some constant /3; typical values of (3 are 2 or 3. Some 
recent work[6l[9] uses a more general power function with minimum regularity conditions, like non- 
negativity, but in all the cases the power function does not depend on the job. Furthermore, most 
energy aware algorithms assume cost functions that are closely related to energy consumption; 
however, in practice the actual energy cost is not simply a function of energy consumption, it 
is a complicated function of discounts, pricing, time of consumption, etc. We consider a more 
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general class of cost functions that are only restricted to be non-negative and can be different for 
different jobs. Because we allow job-dependent energy costs, our algorithms can be used outside 
the CPU-energy setting, where energy cost generally are job independent, and can be applied to 
more general problems that have additional speed-associated costs. Examples of these costs are 
maintenance costs, wear and tear of parts, failure rates, etc. all of which not only depend on 
the speed at which the machine runs, but also the job being processed. We are not aware of any 
other work that allows such general costs. For the weighted tardiness case we require an additional 
regularity condition on the energy cost functions that allows us to control its rate of growth. 

Our paper contains several results for different scheduling problems, we state here the most 
general results: 

Theorem 1.1. Given n jobs with precedence constraints and release dates and a general non- 
negative energy cost function, there is an 0{1)- approximation algorithm for the problem of non- 
preemptively minimizing a convex combination of weighted completion time and energy cost. 

Theorem 1.2. Given n jobs with precedence constraints and deadlines and a general non-negative 
energy cost function, there is an 0{1)- approximation algorithm for the problem of non-preemptively 
minimizing a convex combination of weighted tardiness and energy cost. 

The constants in the 0(1) are modest. Consider the case where we are given a set of speeds 
S = {ai, . . . , am}, at which the machine can run, with aj < (1 + S)aj-i, and some e > 0. Then the 
algorithm for the weighted completion time setting has a 4(1 + e)(l + (5)-approximation ratio when 
only precedence constraints exist, and (3 + 2-v/2)(l + + 5)-approximation ratio when release 
dates are added. The algorithm for the weighted tardiness setting has a 4'^(1 + e)^~^{l + 6)^~^- 
approximation ratio even with arbitrary precedence constraints, where /3 controls the growth of the 
energy cost function. 

1.2 Our Methodology 

The problem of minimizing weighted completion time in the combinatorial setting has been well- 
studied. The work of Phillips, Stein, and Wein [20] and Hall, Schulz, Shmoys, and Wein [171 [IB] 
introduced the idea of a-points, and these have been used in much of the subsequent work. The 
idea is that one first formulates a time-indexed integer program in which decision variable xu is 1 
if job i completes at time t, and then solves its linear programming relaxation. From the solution 
to the relaxation, one computes the a-point of each job, that is, the earliest time at which an a 
fraction of the job has completed in the relaxation. The exact interpretation of when an a fraction 
completes depends upon the particular problem. One uses these a-points to infer an order on the 
jobs and then runs the jobs non-preemptively, respecting that order. There are many variants and 
extensions of these technique including choosing a randomly [131 [T3] or choosing a different a for 
each job [15]. This technique has led to small constant factor approximation algorithms for many 
weighted completion time scheduling problems |24j . 

The time- indexed integer program (IP) formulations for this problem are not typically of poly- 
nomial size. However, the interval-indexed IP, introduced in [18], in which time is divided into 
geometrically increasing intervals and jobs are assigned to intervals rather than individual time 
slots, is of polynomial size. By using this linear program one obtains a polynomial sized linear 
program from which it is still possible to apply the ideas of a-points while suffering only a small 
additional degradation of the approximation ratio. 

In this paper, we extend the interval-indexed IP to handle speed scaling and then design new 
a-point based rounding algorithms to obtain the resulting schedules. In doing so we introduce the 
new concept of a-speeds. We assume, in Sections [21 [3l and [H that we have a discrete set of m 
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allowable speeds S = {(Tj}, and that the rate of power consumption is a polynomial function of the 
speed. In Section 15.31 we describe how to remove these assumptions. Although the time-indexed 
IPs are easier to explain, due to limited space, we will describe only the interval-indexed linear 
programs in this paper. In our interval-indexed IP, a variable Xijt is 1 if job i runs at speed aj and 
completes in interval t. We can then extend the standard interval-indexed integer programming 
formulation to take the extra dimension of speed into account (see Section [2] for details). Once 
we have solved its linear program (LP) relaxation, we need to now determine both an a-point and 
a-speed. The key insight is that by "summarizing" each dimension appropriately, we are able to 
make the correct choice for the other dimension. At a high level, we first choose the a-point by 
"collapsing" all pieces of a job that complete in the LP in interval t (these pieces have different 
speeds), being especially careful with the last interval, where we may have to choose only some of 
the speeds. We then use only the pieces of the job that complete before the a-point to choose the 
speed, where the speed is chosen by collapsing the time dimension and then interpreting the result 
as a probability mass function (pmf), where the probability that the job is run at speed aj depends 
on the total amount of processing done at that speed. We then define the concept of a-speeds, 
which is related to the expected value under this pmf, and run the job at this speed (see Section [3] 
for more details). We combine this new rounding method with extensions of the more traditional 
methods for dealing with precedence constraints and release dates to obtain our algorithms. 

For weighted tardiness, we emphasize again that not much is known about approximating this 
problem, even in the absence of energy concerns. For this problem, we are able to use the same 
interval-indexed linear program, with the objective function modified to tardiness. Because the 
linear program is interval indexed, the non-linear objective function is not a problem. After the 
solving the linear program, we are able to show that with only a constant factor increase in energy 
(over the lower bound from the linear program), we obtain only a constant factor (over the linear 
program) increase in tardiness. Implicit in this analysis is the fact that jobs that receive tardiness 
in the linear program will receive tardiness in our solution; in some sense the speed scaling makes 
accomplishing this easier than in the combinatorial setting. We note that our weighted tardiness 
algorithms does not work in the presence of release dates, as release dates may stop us from being 
able to keep jobs with tardiness in the LP at tardiness in the schedule. 

Finally, in Section \5\ we show how to extend our results for the weighted completion and 
weighted tardiness scheduling metrics to general energy cost functions. We also show how to 
extend our results to the setting where continuous speeds are used and not just a discrete set S, 
while maintaining the same approximation ratio. 

2 Problem Formulation 
2.1 Problem Setting 

We are given n jobs, where job i has a processing requirement of pi E N+ machine cycles, release 
time Tj, and an associated positive weight Wi. Let denote the speed at which job i runs on the 
machine and Cj denote its completion time. Let 11 = {'/r(l), . . . ,7r(n)} denote the order in which 
the jobs are processed, i.e. Tr{k) = i implies that job i is the A:-th job to be processed. Then 
C^(i) = ma'x.{r^n\,C^(i_i)} + is the completion time of the i-th job to be processed, with 
Ctt{o) = 0- We do not allow preemption. 

Let S = {(Ti, . . . ,crm}, be the set of possible speeds at which the machine can run. We will 
assume that crj+i < (1 + S)aj, for some 6 > 0. This is a natural assumption because actual speed 
scaling achieved in CPUs is done via frequency multipliers or dividers. Although a discrete set of 
speeds is probably the most common case for CPUs, in Section [5.31 we show that our algorithm has 
the same approximation ratio when a continuous set of speeds is used. 
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Let Ei{si) denote the energy cost of running job i at speed Sj. For simphcity we initially consider 
Ei{si) = ViPiS^~^, where /3 > 2 and Vi are known constants. As indicated earlier, an energy cost 
function of this form is the standard model for these problems, although our model is more general 
because the energy cost function is job-dependent. In Section [5] we show that our algorithms also 
work for a much larger class of job-dependent energy cost functions. 

The objective is to compute a feasible schedule (11, C), consisting of an order 11 and completion 
times C, possibly subject to precedence and/or release date constraints, and the vector of job 
speeds s = {si, . . . , s„} S WL that minimizes the total cost, 



n 



/(n, S) = ^ [v^PiS^"' + W;.«C,(.)J , (2.1) 
i=l 

Since this function is convex we can assume, w.l.o.g., that each job runs at a constant speed. 

For convenience we will use an extended version of the notation of Graham et al. |16] to refer 
to the different energy aware scheduling problems, i.e. l\ri,prec\ ^ Ei{si) + ^WiCi, will refer to 
the problem setting with 1 machine, with Vi release dates, precedence constraints, and the weighted 
completion time as the scheduling performance metric. Similarly, the l|rj,prec| ^ Ei{si) + ^ WiTi 
will refer to the same setting, but with tardiness as the scheduling performance metric. In all of 
them Ei{si) indicates that the energy cost is also added as a performance metric. 

2.2 Interval-Indexed Formulation 

We now modify and extend the interval-indexed formulation proposed by Hall et al. [18] to accom- 
modate speeds and energy cost. 

The interval-indexed formulation divides the time horizon into geometrically increasing inter- 
vals, and the completion time of each job is assigned to one of these intervals. Since the completion 
times are not associated to a specific time, the completion times are not precisely known but are 
lower bounded. By controlling the growth of each interval one can obtain a sufficiently tight bound. 

The problem formulation is as follows. We divide the time horizon into the following geometri- 
cally increasing intervals: [k, k], (k, (1 + e)/?], ((1 + e)K, (1 + e)^K], . . ., where e > is an arbitrary 
small constant, and k = -^^^ denotes the smallest interval size that will hold at least one whole job. 

' 0"max 

We define interval It = (Tt_i, Tt], with tq = k and Tt = k(1 + e)*"'^. The interval index ranges over 
{!,. . . ,r}, with T = min{[t] : ^(l + e)*-^ > max^^^ + ^lUi and thus, we have a polynomial 
number of indices t. 
Let 

^ _ f li if job i runs at a speed Cj and completes in the time interval It = [rt-i^Tt] ,^ 
^■'^ y 0, otherwise 

By using the lower bounds rj_i of each time interval /j, a lower bound to (j2.ip is written as, 

n rn T 

min^ ^ X] i^iPi^'j'^ + '^i'^t-ij Xijt. (2.3) 
i=l j=l t=l 

The following are the constraints required for the l\ri,prec\ ^Ei{si) + ^WiCi problem: 
1. Each job must finish in a unique time interval and speed; therefore for i = {l,...,n}: 

m T 

j=l t=l 
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2. Since only one job can be processed at any given time, the total processing time of jobs up 
to time interval It must be at most Tt units. Thus, for t = {1, . . . , T}: 



n m t 



i=l j=l u=l ^ 

3. Job i running at speed aj requires ^ time units to be processed, and considering that its 
release time is r^, then for i = {1, . . . , n}, j = {1, . . . , m}, and t = {1, . . . , T}: 

Xijt = 0, ifTt<ri + ^. (2.6) 

4. For i = {1, . . . , n} and t = {1, . . . , T}: 

xu£{0,l}. (2.7) 



5. The precedence constraint ii -< 12 implies that job 12 cannot finish in an interval earlier than 
ii. Therefore for every ii -< 12 constraint we have that for t = {1, . . . , T}: 

m t m. t 

j=l u=l j=l u=l 

It is important to note that this integer program only provides a lower bound for ()2.ip : in fact 
its optimal solution may not be schedulable, since constraints (j2.5p do not imply that only one job 
can be processed at a single time, they only bound the total amount of work in Utit- 

3 Approximation Algorithm for Weighted Completion Time 

We now describe the approximation algorithm for the weighted completion time, called Schedule 
BY a-iNTERVALS AND a-SPEEDS (SAIAS) which is displayed in Figure [3Tl 

Let Xijt denote the optimal solution of the linear relaxation of the integer program (|2.3|) - (|2.8p . 
in which we change constraints ()2.7p for Xijt > 0. In step 1 of the algorithm we compute the optimal 
solution X and in step 2, given < a < 1, we compute the a-interval of job i, which is defined as, 

{m T I 
T :^^Xiju> a\ . (3.1) 
j=l u=l J 

Schedule by a-iNTERVALS and o-speeds (SAIAS) 

Inputs: set of jobs, q € (0, 1), e > 0, set of speeds S = {fii, . . . , cTm}- 

1 Compute an optimal solution x to the linear relaxation (I2.3p - (]2.8p . 

2 Compute the a-intervals and the sets Jj. 

3 Compute an order II" that has the sets Jt ordered in non-decreasing values of t and the 
jobs within each set in a manner consistent with the precedence constraints. 

4 Compute the a-speeds s" 

5 Round down each to the nearest speed in S and run job i at this rounded speed, 

6 Set the i-th job to start at time max{r^(j), where C'"^.^-!^^ is the completion 
time of the previous job using the rounded a-speeds, and C'^^g-j = 0. 

7 return speeds s" and schedule (11°, C°). 

Figure 3.1: Schedule by a-intervals and a-speeds 
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Since several jobs may finish in the same interval, let Jf denote the set of jobs that finish in 
interval It, Jt = {i : = t}, and we use these sets to determine the order 11'^ as described in step 
3. 

Next, in step 4, we compute the a-speeds as follows. Since Yl]Li ^u=i ^iju > ct, we define 
auxiliary variable {xijt} as: 



max |min jxi^vf , a - ^jZl xurf - ft } , o| , t = rf , (3.2) 
I 0, t > rf 



where ft = Ylf=i Ylu=i ^iju < ol. Note that with this auxiliary variable Y^=\ Sm=i — ^- This 
is a key step that allows us to truncate the fractional solution so that for every job i, the sum of Xijt 
up to time interval rf for each speed j can be interpreted as a probability mass function. We define 
this probability mass function (pmf) /ij = (^Uji, . . . , /Ujm) on the set of speeds S = {cJi, . . . , dm} as 



a 



1 

7" ^ '^i-ju- (3.3) 



Let Si define a random variable distributed according to the pmf /Xj, i.e. /Xjj — 
Then, the a-speed of job i, sf , is defined as follows: 



Si 



(3.4) 



We define the a-speeds using the reciprocal of the speeds since the completion times are pro- 
portional to the reciprocals instead of the speeds, and we need to bound completion times in the 
analysis of the algorithm. 

Next, in step 5, because the a-speeds sf do not necessarily belong to the set of possible speeds 
S we round them down to sf, which is the nearest speed in the set such that sf < sf. The following 
lemma bounds the error introduced by this rounding. 

Lemma 3.1. The cost of the solution with the rounded down speeds s"^ is at most (1 + 5) times the 
cost of the solution using the a-speeds s°. 

Proof. The energy cost function Ei{si) is increasing so rounding down does not increase the energy 
cost, but the completion time is now larger. Let be the completion time of job i when the 
speeds s"^ are used and Cf when the rounded ones are used. Since the speeds are reduced at 
most by {1 + 5), then (1 + 5)si" > sf, and we have that, 

Cf = max{r,, (7f_ J + ^ < (1 + J) (^max{r„ Cf„ J + -^^ = {I + 5)Ct , (3.5) 

which implies that Yll=i "^i^f < (1 + <5) Yll=i '^i^f and proves the lemma. □ 

Finally, in steps 6 and 7 we compute the completion times given the calculated speeds and 
return the set of speeds and the schedule (H", C"). 

We now analyse this algorithm's performance for different energy aware scheduling problems. 
In the following subsections we will assume w.l.o.g. that < r^* < . . . r". 
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3.1 Single Machine Problem with Precedence Constraints 

We first need to prove that the output of the SAIAS algorithm is indeed feasible. 

Lemma 3.2. If ii -< i2, then constraint Ii2.8\) implies that < t^. 

Proof. Evaluating the LP constraint ()2.8p corresponding to ii -< i2, for t = t^, we have that, 



j = l 11=1 j = l u=l 

where the last inequality follows from the definition of r^. The chain of inequalities implies that 

ET=l EJil ^^lju > a, so < T^^. □ 

Since the SAIAS algorithm schedules jobs by first ordering the sets Jt in increasing order of t, and 
then orders the jobs within each set in a way that is consistent with the precedence constraints, by 
Lemma [3.2l it follows that the SAIAS algorithm preserves the precedence constraints, and, therefore, 
the output of the algorithm is feasible. Next, we can prove the following result. 

Theorem 3.1. The SAIAS algorithm with a = is a 4(1 + e)(l + S)- approximation algorithm for 
the l\prec\ Y^Ei{si) + ^WiCi problem, with Ei{si) = ViPiS^~^ . 

Proof. Let x*^^ denote an optimal solution to the integer problem (I2.3p - (|2.8|) . xijt the fractional 
solution of its linear relaxation, and xiju the auxiliary variables calculated for the SAIAS algorithm. 
Since in (12. 3p the completion time for jobs completed in interval It is r^-i, it follows that. 



n m T 



n m T 



i=l j=l t=l i=l j=l t=l 



2=1 



The energy terms of the algorithm's solution are bounded as follows. 



ax/3-l 



ViPi 



< ViPiE 



-(/3~1) 





r 11 




1, 
I' 




_Si_ 



-{(3-1) 



ViPiE 



-/3-1 



m 

ViPi^PijO' 



/3-1 
j ' 



(3.6) 



(3.7) 



where the inequality follows from Jensen's Inequality applied to the convex function j/j^. Using 
the definition of pij in (j3.3p and given that 0<a<l,e>0, and Xijt < Xijt, it follows that. 



a{l — a) 



ViP, 



m T 

EE 



13-1- 



j=l u=l ^ ' j=l u=l 

Since there are no release date constraints there is no idle time between jobs. 



(3., 



1 



1 



1 



n m 'i 



j=l 1=1 u=l ' 



j=l J j=l 

and from constraint (12. 5p for t = rf we get, Cf < ^T^rf • 



J = l 1 = 1 u = l ' 



(3.9) 
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:i Tt-iXijt denote the optimal fractional completion time given by the optimal 

solution of the relaxed linear program (j2.3p - (j2.6p . Since it is possible that EJLi Et=i ^ijt > '^'i 

define ^ = a - EJli E[=i"^ ^ijt and ^ = Y.]Li E[=i ^ijt - thus ^ + xf ^ = J2f=i Xijr^ , 
and we can rewrite 

m T-^-l m T 

Ci = Y.Y. ^t-^^nt + ^rf-l^?^ + Trf-xXf^ + 5Z 2^ Tt_lX,,t, (3.10) 
j=\ t=l j = lt=rf + l 

and eliminating the lower terms of the previous sum we get that, 

m T m T 

Ci>Trf-\Xf^^^Y^ ^ rj_iXijt > Tr^-iXf ^ + ^ ^ r^^_i%t = rT-|^_i(l - q). (3.11) 

i = lt=Tf+l j = lt = Tf+l 

Because Trf = (1 + e)Trf -i, from and (fSHT) we get that Cf < ^^^Cj ^ Ya=\ '^i'^f < 
Er=i this, dlS]) and dHS]) it follows that, 

n m T n 



i=i i=i ^ ' 



i=i j=i t=i i=i 



(3.12) 



and we set a = argmino<ct<i | a{i-a) } ~ I' ^° minimize the bound. By Lemma [3. H which bounds 
the final rounding error, we get the desired approximation ratio. □ 

3.2 Single Machine Problem with Precedence and Release Date Constraints 

We now analyse the case with precedence constraints and release dates. Release dates makes the 
problem somewhat harder since they can introduce idle times between jobs. 

Theorem 3.2. The SAIAS algorithm with a = \/2 - 1 is a (3 + 2-v/2)(l + e)(l + 5) -approximation 
algorithm for the l\ri,prec\ ^Ei{si) + 'Y^WiCi problem, with Ei{si) = ViPiS^~^. 

Proof. The bound for the energy terms computed in equation (j3.7p are still valid when there is idle 
time between jobs, we have that, 

, a\B-l ^ (1 + e) V^V^ 13-1- + V^V^ ^QiQ\ 

""'P'^'^^ - ail-af '^'YY''^ ^-^-"^ ail -a) ^^^^Z^Z^^i ^^'^^^ 

j=l u=l j=l u=l 

When bounding the completion time , given the sorting done in step 3 of the SAIAS algo- 
rithm, now one has to consider all the jobs up to the ones in set J^-", and thus, 

< max r,- + V ^. (3.14) 

-ie{Ji,...,j,.} ' , s° 

I 

Since all jobs that have been at least partially processed up to time interval It need to be 
released before rt, it follows that maXjg|jj j^^j rj < Tr^. On the other hand, we also have that, 

i 

t" 1-°' 

E 14 E EE^%^.£iEEE^**.si.,. (3.15) 

je{Ji,...,J^c} 3 j£{J^,...,J^c.} l=lu=\ j=ll=lu=l 
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where the last inequality follows from constraint (|2.5p with t = r". Thus, < 
Ci = X^Jli Ylt=i 'Tt-iXijt, (|3.1ip is still valid and because t^-q = (1 + e)rT-a_i, we get, 

ail — a) ^ a(l — a) ^ 

1=1 

Finally, from (IXTall and (f3J6]l it follows that. 



. Since 



(3.16) 



i=\ 



n m T n 



i=i j=i t=i 



i=l 



(3.17) 



and by setting a = argmino<o<i | ^(1*1^) | — — 1) and again using Lemma l3.ll to bound the 
speed-rounding error, we get the required approximation ratio. □ 

If no precedence constraints and release dates exist, there are two versions of this problem that 
can be optimally solved in polynomial time: when all weights Wi are equal, and when all jobs are 
of the same size (i.e. pi = p, Vi) and all jobs have the same energy cost function. For these cases 
we have the following result: 

i_ 

Theorem 3.3. If Wi = w, Vi or pivf = Vi then the order U is optimal if 



w. 



7r(i) 



> 



-, Vi G {1, . . . , n — 1}. 



Proof. For simplicity we will define S,i = piV^ , q = and Ki = 



(fc-i)- 



^TTT- First, dual formulation 



of problem (j2.ip with no precedence or release date constraints is given by, 



minF(7r) = min^3C^^(j) 



U) 



(3.18) 



i=l 



We now prove both cases by contradiction using the dual formulation. 

When Wi = w, Vi, Theorem 13.31 implies that in the optimal order i^7r(«+i) ^ ^nii)- By contradic- 
tion, let vr be an optimal order such that for some index k, '^7r(fc+i) < '^7r(fc)- For this order the total 
cost is 



{ S.7r{k)in - k + ly + £,^(^k+i)in - ky + Yl C7r(i)(n - i + 1) 

Let TTfc define the order where we switch jobs k and k + 1 from order vr, i.e. TTk{k) = 7r{k + 1) 
and TTk{k -f 1) = TT{k). Given this order we have that 

F{7t) - F{7Tk) = {e.(fc)(n -k+lY + e,(fc+i)(n - kf - e.(fc+i)(n - k + If - ^,(,)(n - kf] 

= ■Kw'i {(n - /c + l)'?(C^(fc) - e7r(fc+i)) -in- k^i^^^^) - ^^(fc+i))} , 
= { {L(k) - L(k+i)) {{n-k + ly -in- ky)} . 
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By our initial assumption the first term is positive (since < S,TT{k)) ^^^d the second one is 

always positive, hence -F(vr) — F{'Kk) > which is a contradiction, since that implies that ttj. has a 
smaller cost. 

For the case when = ^, Vi, Theorem 13.31 implies that an order vr is optimal then i«7r(i) ^ 
^7r(i+i)- Let TT be an optimal order such that for some index k, w^i^i^) < u^7r(fe4-i)- The total cost for 
this solution is 

n \'' ( k / n \^/" \^ n / n 

j=l \i=i / [i=l \j=i J \j=k+l J i=k+2 \j=i 

{k / n \ / \^ n / n 

Y Y'^-ii) + ^-(fc+i) + "^-« + Z]^-« 
i=i yj=j y y j=k+2 J i=k+2 \j=i 

Let TTfc define the order where we switch jobs k and k + \ from order vr. Given this new order 
we have 

F(7r) - F{'Kk) = "I w^f^kj^i) + ^ I - w^(fc) + ^ j > > 0, 

[ \ i=k+2 ) \ j=k+2 J J 

since ^^^^(fc+i) > 'W^7r(fc) by our initial assumption, which is a contradiction since this result implies 
that order vr^ has a lower cost. □ 

4 Extension to the Weighted Tardiness Problem 

In this section we extend our results to the weighted tardiness setting. We still allow for arbitrary 
precedence constraints but no release dates. In this case, each job i also has a deadline dj. The 
tardiness Tj of job i is defined as Tj = max{0, Cj — di}, and the objective function is now given by, 

n n 

^(n, = Y '"iPi^i"^ + Y ^""(i) (^'^c*) ~ ^'^c*)) ■ ^^-^^ 
1=1 1=1 

We now formulate the problem using a modification of the interval-and-speed-indexed formula- 
tion presented in Section [2j Because the completion time can be bounded by X^jLi Ylt=i Tt-iXijt, 
we can bound (|4.ip from below by the following optimization problem, 

n m T 
1=1 j=l t=l 

together with constraints ()2.4p - (j2.8p from the interval- indexed formulation. Note that although the 
objective (j4.ip is non-linear, because we have a interval-indexed formulation, (|4.2p is linear. 

We approximately solve (j4.1j) using the Schedule by o-intervals and q-speeds for Tardi- 
ness (SAIAS-T) Algorithm displayed in Figure HTTJ The main difference with the SAIAS algorithm, 
is that in step 4 we scale up the a-speeds. This scaling makes the completion time of the relaxed LP 
comparable to the completion time of the algorithm's output, and thus jobs that have tardiness 
in the LP also have tardiness in our algorithm. If we rounded speeds down, jobs with tardiness 
in the LP could, at a lower speed, miss their deadline, and thus the approximation ratio could be 
arbitrary large. 

We now analyse the algorithm assuming w.l.o.g. that < < . . . < Since Lemma 
13.21 remains valid, arguments identical to those in Section [2] show that the output of the SAIAS- 
T algorithm is feasible; thus, we have the following theorem: 
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Schedule by a-iNTERVALS and q-speeds for Tardiness (SAIAS-T) 

Inputs: set of jobs, q € (0, 1), e > 0, 7 > 1, set of speeds S = {ui, . . . , am}- 

1 Compute an optimal solution x to the linear relaxation (I4.2p . ()2.4p - ()2.8p . 

2 Compute the a-intervals and the sets Jt as in the SAIAS algorithm. 

3 Compute an order 11" that has the sets Jt ordered in non-decreasing values of t and the 
jobs within each set in a manner consistent with the precedence constraints. 

4 Compute the a-speeds s" and scale each sf to sf = jsf. 

5 Round up each sf to the next speed in S, sf and run each job i at this new speed. 

6 Set the i-th job to start at time max{r^(-j), C'"^._-^j}, where is the completion 
time of the previous job using the rounded a-speeds, and ^^"(o) ~ ^• 

7 return speeds s" and schedule (H", C"). 



Figure 4.1: Schedule by a-intervals and a-speeds for Tardiness Algorithm 
Theorem 4.1. The SAIAS-T algorithm with 7 = ^^±^ and a = \ is a 4^(1 + ef-^{l + Sf'^- 
approximation algorithm for the l\prec\ -|- Y^WiTi problem, with Ei{si) = ViPiS^~^ . 

Proof. Let Ci = X^JLi Ylt=i Tt-iXijt denote the optimal fractional completion time of the relaxed 
linear program. {Ci — di)~^ is a lower bound for the optimal tardiness (C^* — di)~^, since ^jt{Tt-i — 
di)^Xijt > {Ci - di)'^. Thus, 

n m T n n m T n 

^^Yl '"iPi'^j'^^ijt + Y'^i {Ci - diy <YYY^ ViPia^~'^xy + ^Wi{Ci - di)^ . (4.3) 
i=l j=i t=l i=l i=l j=l t=l i=l 

Let denote the completion time of job i using speeds s" and C° the one using speeds s". 
Because there are no release date constraints, there is no idle time in between jobs; therefore, 

^r = t§^='-t^='-cr. (4.4) 

j=l ^ j=l '3 T 

Since ^Mi remains vahd, it follows that Cf < ^ < ^^^Q- The key step is 



a) ^ i — 7 ol{1—ol) 

a(l-o) ' 



that by setting 7 = , which makes the two completion times comparable, we have that. 



5^u;.(cr-<i.) <E^^[-^7^^C.-dA =Y,w.{C.-df)\ (4.5) 

i=l i=l \ I \ / / i—i 



The energy term is bounded in a manner analogous to 

m T 



ViP^{sff-' = i^-'vip.{sff-' < il±il!^^,p,f; J^af-^x.,,, (4.6) 

(a(l - a)Y t=i 



where the last inequality follows from (j3.8p that remains valid. 
From ([O]), dlS]), and ([13]) it follows that. 



i=i i=i ^ ^ {a{l-a)Y 



n m T n 
i=l j=l t=l i=l 



11 



Energy Aware Scheduling v.arXiv 4.6 



Carrasco, Iyengar, Stein 



Because speeds are rounded up, the completion times, and thus the tardiness can only improve, 
whereas the energy cost increases. Since at most we speed up each job by a factor (1 + 6), we have 
that, 

n n 

E,{st) < E, ((1 + 6)s't) = (1 + 5f-'E,{sf) Ei{-sf) < (1 + 5f~' ^^(^°)- (^.7) 

1=1 i=l 

The approximation ratio follows from setting a = argmino<Q-<i | ^^^^^ ^ a))^ } ~ ^' ^^^^"^^y 

could use J^(^'^^-^ti-^ (|4.6|) to compute a tighter bound, but the resulting expression is not as 
simple. □ 

We are not able to extend this algorithm for the V\ri\ ^Ei{si) + ^WiTi problem, since it is 
based on speed scaling to make sure that jobs are finished within a desired time interval. When 
release dates are present, we do not see how to arbitrarily reduce the completion times. 

5 Extension to General Energy Cost Functions 

In this section we consider the extension to general energy cost functions, as opposed to simply 
energy consumption. We begin by considering discrete speeds, as in the previous sections, but in 
Section [5.31 we will relax this requirement. 

Managers of data centres are clearly interested in the energy cost metric, since they need to 
balance the penalty for violating the service level agreements with the cost of energy. The energy 
price curves for industrial consumers are often quite complicated because of energy contracts, 
discounts, real time pricing etc.; therefore it is very important to consider general cost functions 
in the scheduling model. Hence, in this section we use £i{si) as the general energy cost function 
of running job i at speed Sj. We will require that £i{si) is non-negative, just as in O [9], but no 
other requirements are needed for the weighted completion time setting. For the weighted tardiness 
setting we will require an additional regularity condition that bounds the growth of the energy cost 
function. 

Since in practice the processor speed can be dynamically changed during the course of a job, 
one can replace the general cost function by its lower convex envelope. Hence, without loss of 
generality, we can assume that £i{si) is convex. Furthermore, since the machine can only run at 
the speeds in S, we can also consider that £i{s) is linear in between these speeds. Hence, for every 
s S [(Tj,(Tj+i] such that s = Acjj + (1 — A)(Tj+i, with A G [0, 1], then £i{s) = \£i{a j) + {1 — X)£i{a j+i) . 

Note that for bounding the energy cost terms in the weighted completion time setting, we only 
used the fact that the energy consumption function Ei{s) = ViPiS^~^ is convex. Thus, the previous 
bounds extend to our more general class of functions £i{s). In the weighted tardiness case we 
required also a bound on the growth of the energy cost function, which we will address in Section 
K2[ 

5.1 Weighted Completion Time Problem with General Energy Cost 

The objective function (|2.3p is extended as follows, 

n m T 

min^^^ (£:j(f7j) + WiTt-i) Xijt, (5.1) 
i=l j=l t=l 

where £i{crj) are just coefficients. Given that we only change the energy cost related terms, all the 
completion time related bounds computed previously are still valid. 



12 



Energy Aware Scheduling v.arXiv 4.6 



Carrasco, Iyengar, Stein 



The only modification required is in the rounding procedure at the end of the SAIAS algorithm, 
where it was done by rounding down the a-speeds. Now instead we will round them up or down 
such that £i{sf) < £i{sf), which is always possible since £i{si) is linear in between the speeds in S. 
With this change Lemma l3 . 1 1 remains valid and we can extend the algorithm to our general energy 
cost functions. 

Theorem 5.1. The SAIAS algorithm with a = ^ is a 4(1 + e)(l + 5) -approximation algorithm for 
the l\prec\ ^£i{si) + ^wiCi problem, for all general non-negative energy cost functions £i{s). 

Proof. Because £i{a), i = {l,...,n} are convex functions, (|3.7p remains valid since £i{sf) = 
£i(E[si]) < ¥,[£i{si)] = Aiij'£'i(o"j), and thus, from the definition of fuj, and from < a < 1, 

e > 0, and Xijt < Xyt, 



n m T 



< ^EEE^^(^.)^^.-* ^ i^EEE^^(^^)%*- (5.2) 

i=i i=i j=i t=i ^ ' 1=1 j=i t=i 

The proof follows since the bounds for the completion time in Theorem 13.11 remain valid, as 
well as Lemma |3. 11 □ 

By the same argument we also have that. 

Theorem 5.2. The SAIAS algorithm with a = \/2 — 1 is a (3 + 2^/2)(l + e)(l + 5) -approximation 
algorithm for the l|rj,prec| ^ ifj(sj) + ^WiCi problem, for all general non-negative energy cost 
functions £i{s). 

5.2 Weighted Tardiness Problem with General Energy Cost 

We replace the energy term in (14. 2p with the general energy cost term to obtain the new objective 

n m T 

minEEE + ^« ("^i-i ~ '^^)^) ^vt- (5-3) 

i=l j=l t=l 

Since the SAIAS-T algorithm speeds up the jobs, we need to add the following regularity 
condition for the energy cost functions £i{cr) in order to obtain performance bounds: 

Assumption 5.1. 3/3 G N"*", such that £i{'yai) < ^l^~^£i{ai), V7 > 1. 

Theorem 5.3. The SAIAS-T algorithm with 7 = -^^^ and a = \, is a^P{l + ef-^{l + Sf'^- 
approximation algorithm for the l\prec\ £i{si) -\-^WiTi problem, for all non-negative energy cost 
functions £i{s) that satisfy Assumption \5. 1[ 

Proof. As before, all the completion time related bounds (14. 4p and ()4.5p remain valid, so only a 
bound analogous to (|4.6p is needed. From Assumption 15.11 it follows that, 

£^{sf) < iP-%{sf) < A-t£i^ J^£:,(^^.)%*- (5.4) 

^ ' i=i t=i 

Thus, from (j4.5p it follows that, 

n m. T n 



i=l 1=1 



E E E ^^(^iXt + E - 
i=i j=i t=i 1=1 



Since we are rounding speeds up, equation (|4.7p remains valid and thus taking a = ^ completes 
the proof. □ 



13 



Energy Aware Scheduling v.arXiv 4.6 



Carrasco, Iyengar, Stein 



5.3 Continuous Speeds 

As commented previously, our algorithms are also applicable for the case when a continuous set 
of speeds is possible. In this case we modify the SAIAS and SAIAS-T algorithms, eliminating the 
rounding step required at the end of each algorithm. 

When the operating range of the machine is given, i.e. the speed limits cJi„in and (Tmax; since 
our IP requires a speed index, we need to quantize the set [cmim c"max] in different speeds. We 
can do this by setting ai = cTmin, and as before we define speed aj = (1 + 6)aj-i, for some 6 > 0, 
making sure that am > o"max in order to cover the whole operating range. Just by rounding as 
described in Section [5. II for the weighted completion time setting and rounding up for the weighted 
tardiness setting we can prove the following lemma: 

Lemma 5.1. The optimal solution for the IP i2.3\) - [278\) is at most (1 + 6) times the optimal 
solution of the energy aware problem in the weighted completion time and continuous speed setting, 
and the optimal solution for the IP ^.2^ , Ii2.4\ )- f2^) is at most {l+6)^~^ times the optimal solution 
of the energy aware problem in the weighted tardiness and continuous speed setting. 

The proof is similar to Lemma 13.11 for the weighted completion time and similar to equation 
(j4.7p for the weighted tardiness setting. 

Since there is no additional rounding at the end of the algorithm, using Lemma l5. II we get the 
same approximation ratios as in Theorems 15.1} 15. 2|, and 15.31 

When the operating range of the machine is not given, and we are interested in determining 
a set S that covers the optimal speeds from the continuous case, we need the following additional 
regularity condition on the energy cost functions: 3^ < oo such that £j{si) is increasing Vsj > ^. It 
is easy to prove that this is a necessary and sufficient conditions for the problem to be well defined, 
and thus we can compute iTmin and <Tmax such that the optimal speeds s| € [o"miin <7max] > for all 
i. Then we can apply the same procedure as before to quantize and build the set of speeds, and 
proceed to compute an approximate solution. 

6 Conclusion 

In this work we described new techniques for developing constant approximation algorithms for 
energy aware scheduling problems with very general job-dependent energy cost functions, that 
work on both discrete and continuous speed sets. Furthermore, we present the first algorithm, to 
the best of our knowledge, that tackles the energy aware weighted tardiness setting, even in the 
presence of arbitrary precedence constraints. 

We believe that our methodology, which extends the idea of Q-points to the energy aware setting 
by developing the a-speeds concept, should have many more applications. We suspect that, via 
techniques such as using randomly chosen values of a or using different a values for different jobs, 
we could obtain tighter bounds, and also that these techniques could be extended to other settings, 
such as multiple parallel machines among others. 
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